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Abstract 

In the framework of the quasiparticle-phonon model, we study the non- 
statistical proton decay of excited states in odd nuclei towards low-lying col- 
lective states. Partial cross sections and branching ratios for the proton decay 
of the high angular momentum states in ^^Sc, ^^Cu and ^^Nb are evaluated. 
The calculated branching ratios predict strong direct proton decays to the 
low-lying vibrational states in ^"^Sc and ^^Nb. A general agreement with ex- 
isting experimental data is found. 
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I. INTRODUCTION 



Nucleon transfer reactions induced by hadronic probes at intermediate energies favour the 
excitation of high angular momentum states lying above particle emission threshold 
The observed structures originate from the coupling of the initial single-particle mode with 
more complex states. The coupling of single-particle states with surface vibrations is mainly 
responsible for the damping of the single-particle mode • The particle decay of highly 

excited states gives the opportunity to study in detail the damping process. For example, 
the relative contributions of the direct and statistical components to the damping of single- 
particle mode can be found. Up to now, experimental data and theoretical calculations are 
available mainly for the neutron decay of high-lying single-particle modes. Very recently, 
the proton decays of high-lying states in "^^Sc, ^^Cu and ^^Nb have been measured 

In odd nuclei, the simplest excited states can be described as admixtures of single- 
particle states or quasiparticle states coupled to collective excitations (or phonons) of the 
even-even core. This weak coupling picture has been successfully applied to obtain the 
strength functions of a variety of odd nuclei [|l],^ . A method for calculating particle escape 
widths in the framework of the quasiparticle-phonon model (QPM) has been suggested 

in Ref. where the inclusive semi-direct neutron decays of high-lying states in 209p]^ have 
been studied using the general procedure of Refs. [jrT]-^ which was established for the case 
of nucleon emission from giant resonances. In a previous work we have extended the 
method of Ref. ^ in order to calculate non-statistical particle decays of excited states in odd 
nuclei leading to exclusive channels which correspond to the ground and low-lying excited 
states of the even-even core. In the present work we apply our method to study the proton 
decay. Using the QPM we calculate the partial cross sections and branching ratios for the 
proton decay of the high angular momentum states in ^^Sc, ^^Cu and ^^Nb and compare 
them with experimental data. 

This paper is organized as follows: in Sec. II we describe briefly our theoretical approach 
to treat the direct nucleon decay of high angular momentum states of single-particle type. 
In Sec. Ill a comparison of the calculated and measured branching ratios for non-statistical 
proton decay of high- lying states in "^^Sc, ^^Cu and ^^Nb is presented. Finally, in Sec. IV 
conclusions are drawn. 



II. THEORY 



A. The projection operator method 



The projection operator method is a convenient approach to treat problems involving 
single-particle continua in a many-body context. The general formalism was introduced by 



Yoshida and Adachi |TT[] and it has been applied to studies of high-lying states in odd nuclei 
PJT^. The detailed expressions can be found in Ref. Here, we simply recall the main 
features. 

We write the hamiltonian of the A+1 system in the form: 



H = h + Hcore + H, 



coupl 



(1) 
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The first term describes the motion of a particle in an average potential U created by the 
particles in the core: 

A^-^V^ + C/. (2) 

The core hamiltonian is a sum of single-particle hamiltonians hi and two-body residual 
interactions Vij\ 

A A 

Hcore = X! + X! • (^) 

i=l i<j 

The last term of if is a sum of interactions between the odd particle and the core particles: 

A 

Hcoupl = X! ^o,i ■ (4) 

i=l 

The physical spectrum of h consists of a small number of bound states {ipi,ei} and a 
continuum of scattering states {(pe, e} which form altogether a complete, orthogonal basis. 
The projection operator method consists in introducing another complete set of orthogonal 
basis states which is a direct sum of two complementary subsets, a first subset of discrete, 
orthonormal states {(pa, ^a} which span the single-particle q space and a second subset of 
{(/>£, e} continuum states spanning the complementary p = 1 — q space. We will denote by 
al^,aa {cil,a^) the creation and annihilation operators of state ipa (v^J. 

In the spirit of the QPM ||^ the hamiltonian Hcore is treated in the random-phase ap- 
proximation (RPA) in a discrete space, i.e., the particle- hole configurations of RPA are built 
only with q-space states. We denote by and Ol the energies and creation operators of 
these RPA states which describe core excitations. If | > represents the RPA ground state 
of the core, the properties of the (A-l-l)-nucleus except for its nucleon decay properties can 
be described in terms of the one-particle states al^ \ > and one-particle-plus-phonon states 



a 



i ® ol] I >. We can write: 



di> = d| I > 

EcM + E<>^®ot]) |o> . (5) 



We call Q space the space spanned by the (real) state vectors | di > and Q the corresponding 
projection operator. The amplitudes Cj^^ and , and the energies Ui oi \ di > are 
determined by diagonalizing H in the RPA, i.e., one solves: 

[QHQ,4]=u4, (6) 

within the approximation of commutator linearization. The distribution of | C*^*-* ^ repre- 
sents the strength function from which one can deduce the spectroscopic factors 

To allow for nucleons to decay it is necessary to introduce state vectors where the odd 
particle has a non-zero probability of being at infinity. This is achieved by constructing the 
P space complementary to Q space and consisting of all states which are linear combinations 
of the following one-particle and one-particle-plus-phonon configurations: 
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I e >= a,^ I > , \e,iy>= [al ^ Ol] \ > . (7) 

The present definition of P space neglects continuum effects on tlie plionons Ol which 
can also in principle couple to the continuum and emit nucleons by themselves. Actually, 
these effects should have a small influence on the non-statistical particle decay of the (A+l)- 
nucleus since the most important phonons contributing to the particle-phonon coupling are 
the low-lying collective states of the core. 

The direct sum of Q space and P space is by construction the complete particle-plus- 
phonon space in which the hamiltonian H should be solved. It is completely equivalent to 
solve in the simpler Q space the more complicated effective hamiltonian: 

n{E) ^ QHQ + QHP ^^^^^p^p PHQ 

= Hqq + W{E) , (8) 

where P is the projection operator onto P space {P + Q = 1) and E is the energy of the 
system. The hamiltonian H is complex and energy dependent. For each value of E one has 
to find the set of complex states and eigenenergies: 

I '^^) ^ ^} I 0) , 

= a;, - '-Tj , (9) 

satisfying: 

[n{E)M] = ^^1^^ ■ (10) 



B. Escape widths 

We consider a direct transfer reaction a + A b+{A+l)* followed by a sequential decay 
(^4-1- 1)* — > p + A* where the (A-l-l)-nucleus in a highly excited state decays by a semi-direct 
proton emission. In this process, the | Vi > states will act as doorway states. If we describe 
the reaction mechanism in a simple approach, e.g., a distorted wave Born approximation 
(DWBA), we can write down the scattering amplitude from an initial channel i where the 
target A is in its ground state | > to a final channel / where the residual nucleus is left 
in some excited state | u >= Ol \ > with excitation energy Ei, while the escaping proton 
has an energy E — E^. Using the complex bi-orthogonal basis {| Vi >,<'Di |}, we have: 

y _ y^y^ < <pi;\E-E,),u\ H\V,><Vi,b\V\0,a> 

where V is the interaction inducing the particle transfer from a to A, (p^^^ is an incoming 
wave of p space at energy E — E^, and the sum over is restricted by the angular 

momenta and parities of states | T>i > and u. The case where the final channel is the ground 
state of the residual nucleus corresponds to the above expression with i/ — 0, E,^ — 0. 
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The nucleon-transfer matrix element < T>i,b \ V \ 0, a > is proportional to the one- 
quasiparticle amplitude of the state \ Vi >. In analogy with the amplitude C*-*^ of Eq.(^), 
we denote it by C^'^\ We also introduce the partial escape amplitudes of state \ > to 
channel 

-fi,,{lj) = V2^<4i\E-E,),u\H\Vi> . (12) 



The partial widths are: 

rl = E 1 P • (13) 

ij 

To obtain an expression for the cross section simple enough to lend itself to a discussion in 
terms of escape widths, let us assume furthermore that interference terms between different 
doorway states can be neglected. The density of \Vi) states is large, and for each interval 
centered around E and containing states we define locally averaged quantities: 



rI(E) = Erl/iv 



T\E)^Y.^l{E) . (14) 

V 

Then, one can rewrite the cross sections in the following form: 

a,(E) oc C\E)Vl[E)lV\E) . (15) 
The branching ratios can be calculated by the following formula: 



C. Inputs of the model 

The above formalism is applied to study the semi-direct proton decay of the nuclei "^^Sc, 
^''Cu and ^^Nb with the aim of comparing the predictions with existing data from exclusive 
measurements. 

Although the QPM model Soloviev et al. |?;10] is not fully consistent since the resid- 



ual interaction between quasiparticles is not derived from the quasiparticle mean field, it 
has the advantage that its two-body residual interaction is chosen of a multipole-multipole 
separable form and therefore, it enables one to work with configuration spaces of large di- 
mensions without facing the cumbersome problem of diagonalizing large matrices. This is 
very important for the results to be meaningful as it was already shown in a previous work 
on neutron emission . 



The mean potentials entering the single-particle hamiltonians hi are of Woods-Saxon 



form. Their parameters are chosen according to Refs. ||15|jl6[| with some readjustments 
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such that, once the couphng to phonons is included, the ground states are located at their 
respective experimental positions with respect to the proton separation threshold. The QPM 
hamiltonian includes the monopole pairing which must be taken into account for nuclei with 
open shells This is the case when we calculate the core properties of ^^Cu and ^^Nb. 

The residual particle-hole interaction of Eq. (|^) is taken of a separable form in coordinate 
space with effective interaction strengths considered as adjustable parameters. For the radial 
interaction form factor we have used /(r) = dU/dr where U{r) is the central part of the 
Woods-Saxon potential. For each excitation mode of the even-even core corresponding to a 
given angular momentum, parity and isospin, the interaction strength is found by requiring 
that the lowest collective state calculated in RPA be at the experimental energy. These 
interaction strengths are also used for the quasiparticle-phonon coupling. We have modified 
the single-particle proton spectrum for ^^Nb in comparison with Ref. to get a better 



description of high-lying 1(77/2, l/iii/2 subshells, but so far as energies of levels near the Fermi 
surface are practically the same there are no changes in the properties of low-lying vibrational 



states in comparison with those of Ref. ||I4[. For example, the calculated B(EA) values of the 
first low-lying collective states 3^, 5^ in ^°Ca are B(E3)=1.36e2/"^^ B(E5)=2.67eV^^° that 
can be compared with the corresponding experimental values 1.24e^/m^ and 2.97e^/m^°. 
For all nuclei under consideration the calculated values of B(EA) are in general agreement 
with experiment. In the present model, these low-lying phonons of the core are the physical 
channels where the initial state in the excited odd-A nucleus can decay by semi-direct proton 
emission. In the actual calculation, a very large number of RPA phonons are included in 
the quasiparticle-phonon basis, but for "^^Ca the first quadrupole state cannot be treated as 
a one-phonon or particle-hole state and it is outside the present model. On the other hand, 
the experimental data show that the proton decay of "^^Sc to the 2f state in "^^Ca is very 
weak. 

At high excitation energies, when the level density becomes large, it is often more con- 
venient to calculate the strength function []TB[ instead of solving the secular equations (PUD. 
One defines the strength function C'^{E) as the strength distribution IC^*-*]^ folded with an 
averaging function p : 

cl{E) = Y.\c]^\'p{E-a-;) , (17) 

i 

where p{z) is usually chosen as a Breit-Wigner function, p{z) = ^ ^2^.^2/4 - It is shown 
in Ref. |T^ that the calculation of the strength function (^"^{E) can be done without the 
detailed knowledge of the amplitudes |C^*^|2. The same averaging procedure can be apphed 
to calculate the partial escape widths and the branching ratios. In all calculations performed 
in this work we have adopted the value A = 0.1 MeV. 



III. RESULTS AND DISCUSSION 



Using the methods presented above we have calculated the cross sections and branching 
ratios for the non-statistical proton decay of high-lying states with high angular momenta 
in ^^Sc, ^^Cu and ^^Nb going to the ground and low- lying collective states of ^''Ca, ^^Ni and 
^°Zr, respectively. 
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An example of cross section calculations by the strength function method is shown in 
Fig. 1 for the case of proton decay of the (^9/2 states in ^^Sc. Since the calculated cross 
section is proportional to the square of the one quasiparticle amplitude of the wave function 
this figure enables one to see the lgg/2 strength distribution in ^^Sc, too. The lgQ/2 strength 
is distributed over a broad energy interval due to the coupling with the collective vibrations. 
Results of our calculations for the lgg/2 proton particle strength and experimental data 
obtained from the ^°Ca(^He,dp) reaction at E3He='24:0 MeV []T^ are presented in Table 
1. It is seen from this table that the calculations reproduce reasonably well the integral 
characteristics of the lgg/2 proton strength distribution in ^^Sc. 

The partial contributions of the non-statistical proton decay into the ground state and 
3i, 5i states of ''^Ca are presented in Fig. 2. It is seen from Fig. 2 that the partial cross 
sections are very energy dependent and as a result different channels can dominate at some 
excitation energies. The solid curve in Fig. 1 presents the sum of all 3 partial channels. 
As one can see from Fig. 2 at excitation energies below 8.4 MeV the ground state channel 
dominates, but there are strong transitions to the 5i state at 7.5 and 8.0 MeV. The 3j~ and 
5i channels begin to contribute importantly starting from 8.4 MeV and they become the 
main contributors after 9.5 MeV. 

The sum of partial widths for different channels in some energy intervals are given in 
Table 2. As one can see from Table 2 for the energy interval 2.0 - 12.4 MeV the proton 
decay to the 5j~ state gives a contribution of about 54% and the 3i channel contributes 
about 37% in the total sum. There is a predominance of the 5j~ channel in comparison 
with the 3i one in the energy interval 7.0 - 12.4 MeV. Such a behaviour can be understood 
from the structure of the decaying states and the angular momentum and energy carried 
away by the emitted proton. It is seen from Eq. (0) that the partial escape amplitude 
for fixed values of Ij is proportional to the contribution of the quasiparticle-plus-phonon 
configuration made of a quasiparticle with the same Ij and a phonon of the final state. The 
cross section depends on the square of these amplitudes summed over all Ij allowed by the 
angular momentum coupling rules. In the case of the proton decay of the 5^9/2 states at 
excitation energies 7-12 MeV into the 3j~ channel the outgoing protons carry (/, j) = (l,3/2) 
mainly whereas it is (/, j) = (l,3/2) and (/, j) = (l,l/2) for the 5~ channel. Thus, penetration 
factor arguments favour the latter channel. It is well known that the strongest coupling 
between the single-particle states and phonons takes place for the collective phonons. That 
is why the coupling with the low-lying vibrations is mainly responsible of the damping of 
the high-lying single-particle modes and giant resonances in nuclei. This is the case for 
the 3i and 5j~ states in ^°Ca. By this reason the strengths of the 3J~ and 5j~ channels are 
redistributed in a broad energy interval. The decrease of the spectroscopic strength with 
an increase of the excitation energy leads to a decreasing of cross sections of the channels 
discussed above. It is worth to mention that a very similar behaviour for the 3j~ and 5^" 
channels for the neutron decay of the high angular momentum states in ^°^Pb also takes 



place Our conclusions about the role of different channels for the proton decay in Sc 



agree with observations made in Ref. |T7], but no separation of the direct part from the 
statistical one for the proton decay was done in that work. 

As a second example we consider the proton decay of the gg/2 states in ^^Cu. Results 
of our calculations for the gg/2 strength distribution in ^^Cu and experimental data [jl9| are 



presented in Table 3. The calculations reproduce correctly the strength near 7 MeV and 
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this part is mainly responsible for the proton decay. For the excitation energy interval (2 
- 8.9) MeV in ^^Cu the summed ground state width is equal to 10.3 eV and this is much 
less than in the ^^Sc case. An additional contribution from the 2"*" channel has a summed 
width that is six times less than that of the ground state channel. There are no transitions 
to other vibrational states. This is easily understood if one looks at the structure of the 
states under discussion. For the 2"*" channel the outgoing protons carry (/, j) = (4,9/2) and 
there is a suppression of such transitions due to penetration factor effects. The configurations 
including 3J~ and 5j~ states coupled with continuum are located at somewhat higher energies. 
This is why no proton decay can proceed to the 3j~ and 5j~ states in ^^Ni. 

Very recently, detailed experimental information about proton decay of isobaric analog 
states (IAS) in ^^Nb has been shown in Ref. The IAS are strongly excited by means of 
the ^°Zr(a, t)^^Nb reaction. A sharp peak is seen at 12 MeV excitation energy. Around this 



excitation energy three IAS have been observed earlier ||20[ at 11.93 MeV, 12.07 MeV and 
12.15 MeV. It is shown in Ref. |^ that at 12 MeV by means of {a, t) reaction predominantly 
the /ill/2 state is excited. The proton decays to the ground state and some low- lying states 
are observed. The partial differential cross sections for each final state are given. The data 
reveal proton decays predominantly to 5j~ and 3^ excited states of ^^Zr. 

Experimental data |^ for energies and spectroscopic factors of IAS and results of our 
calculations are given in Table 4. Besides, this table also contains calculated partial widths 
for the proton decay of three IAS to the ground, 2i, 3f , and 5^ states. 

As one can see from Table 4 the calculated excitation energies for the 7/2+ IAS are higher 
than the experimental ones but the spectroscopic factors extracted from experimental data 
within 20% accuracy are reproduced reasonably well. 

Let us discuss the proton partial widths for these IAS. It is seen from Table 4 that for the 
first 7/2+ state the main decay channel is the 4+ channel. For this state the (/, j) = (0,l/2) 
protons are most important for the 4+ channel and there is a rather weak transition for 
the 2f channel due to the (/, j) = (2,3/2) protons. The ground state width is proportional 
to the one quasiparticle strength (the spectroscopic factor). Since the contribution of the 
one quasiparticle component in the norm of the wave function of this IAS is much less than 
contributions of the quasiparticle-plus-phonon components the transition to the ground state 
is weaker than to the 4f. The particularities of the partial proton widths for the second 
7/2+ state can be understood from the structure of this state, too. For the 4^ channel the 
(/, j) = (2,3/2),(2,5/2) and (4,7/2) protons can contribute besides the (/, j) = (0,l/2) protons. 
As a result the decay width to the 4+ state increases strongly for the second 7/2+ IAS 
in comparison with that of the first 7/2+ state. An increase of the contribution of the 
configuration constructed from the first quadrupole phonon and 2(^3/2 single-particle state 
in comparison with their contribution to the structure of the first 7/2+ state results in an 
essential growth of the decay width to the 2+ channel. The ground state width becomes 
larger mainly because of the increased spectroscopic factor. The transitions to the 5^" state 
can take place due to the outgoing protons with (Z, j) = (5,ll/2) but they are supressed 
because of penetration factors and a small contribution of relevant components in the wave 
function structure. 

In the case of the ll/2~ IAS the 5~ channel dominates and this is due to the outgoing 
protons carrying (/, j) = (0,l/2),(2,3/2) and (2,5/2). The components containing the 3^" 
phonon have a small contribution in the wave function structure and as a result the width 
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for the 3j~ channel is small, too. In spite of unfavourable penetration factor the 2i width is 
4 times larger than that of the 3j~ channel. This is due to the configuration consisting of the 
25" phonon coupled with hii/2 which contributes about 15% in the norm of the 11/2^ IAS. 

The comparison of the calculated partial widths with the data of Ref. P] shows that 
there is a qualitative agreement. The main channels of the proton decay of (ll/2)~ IAS are 
reproduced in the calculated structure. The dominance of the 5^ channel is well established 
and the calculated partial width of 3.3 keV (Table 4) is in agreement with the measured one 
(2.9 keV). The 3j~ channel is more pronounced than ground state channel but the calculated 
partial widths for both channels are much less than the measured ones. The calculations 
indicate a large decay to the 2f channel. Such decay is discussed in Ref. but quantitative 
estimations for the contribution of the 2^ channel in the cross section have not been evaluated 
in that work. 

The summed partial widths for the proton decay of the three IAS to the low-lying 
vibrational states are also presented in Table 4. As one can see from the last row of Table 
4 the 4^ channel exhausts almost 60% of the total cross section, the 5j~ channel contributes 
about 26% and the 2f channel gives about 10% of the total strength. The contribution of 
the ground state channel for the proton decay of the three IAS in ^^Nb is 4% only. 



IV. CONCLUSIONS 

A microscopic approach based on the QPM has been applied to calculate the non- 
statistical proton decays of high angular momentum states excited in one-nucleon transfer 
reactions. Partial cross sections and branching ratios for proton emission from high-lying 
states in '^^Sc, ^^Ni and ^^Nb have been evaluated. The calculated branching ratios demon- 
strate the existence of strong direct proton decays to the low- lying vibrational states in ^^Sc 
and ^^Nb and enables one to understand particularities of the decay to different channels. 
One can conclude from an analysis of calculated partial cross sections that, for high angular 
momentum states the non-statistical proton decay is more favourable into the higher angular 
momentum and lower excitation energy final states when the two following conditions are 
fulfilled: a strong particle-phonon coupling in that channel and a penetration factor which 
is not hindered by angular momentum or energy. A similar conclusion was also reached 
for the neutron decay case [|T^. A general agreement with existing experimental data is 



obtained and the predicted branching ratios for different channels can be used to analyze 
future experimental data. 
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FIGURES 

FIG. 1. Cross sections (arbitrary units) for the excitation and proton decay of 39/2 states in 
^^Sc into low-lying states, as a function of excitation energy E. The solid line is the sum of ground 
state, 3^ and 5^ channels, the dotted line shows the partial contribution of the ground state 
channel. 

FIG. 2. Partial cross sections (arbitrary units) for the excitation and proton decay of gg/2 states 
in ^^Sc as a function of excitation energy E. The solid line is the 3^ channel, the dotted-dashed 
line is 5]^ channel and the dotted line is the ground state channel. 
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TABLES 



TABLE I. The (^9/2 proton particle strength in Sc . 



AE^(MeV) 


E^iexp.) (MeV) E^{calc.) (MeV) C^S{exp.) 


C^Sicalc.) 


2.0 - 12.4 


8.80 8.96 0.84 


0.92 


8.4 - 12.4 


10.4 9.49 0.56 


0.76 




5.04 4.65 0.12 


0.05 


TABLE IL Summed partial widths (in MeV) for different channels for the proton 


decay in ^^Sc 


AE^(MeV) 






2.0 - 7.2 


0.225 3.63 10-5 


2.14 10-2 


2.0 - 8.4 


3.314 0.170 


2.31 10-2 


8.4 - 12.4 


1.689 11.36 


7.832 


2.0 - 12.4 


5.003 11.53 


7.855 


TABLE in. The 1^9/2 proton particle strength in ^^Cu . 


^^^(MeV) 


C^S{exp.) 


C^S{calc.) 


3.0 


0.35 


0.60 


~7.3 


0.35 


0.35 


3.1 - 7.2 


0.09 


0.05 



TABLE IV. Partial widths (in keV) for different channels for the proton decay in ^^Nb 



E:, (rj-/).) (MeV) E,, (r-r/Zr.) (McV) 


.V 


C^Sierp.) 


C-S(calc.) 


-pt pT pt pT 
2,1 :-!, 1,1 




11.79 12.00 


7/2+ 


0.11 


0.18 


0.156 0.019 0.990 


0.005 


12.07 12.15 


11/2- 


0.35 


0.57 


0.049 0.296 0.073 0.002 


3.308 


12.15 12.90 


7/2+ 


0.40 


0.29 


0.309 1.010 6.738 


0.027 








Er^ 


0.514 1.325 0.073 7.73 


3.340 








% of total 


4.0 10.2 0.6 59.5 


25.7 
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